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Application  of  extreme  value  statistics 
to  corrosion  is  reviewed.  It  is  empha- 
sized that  the  concept  of  corrosion 
probability  is  important  for  a quantita- 
tive evaluation  of  corrosion  failure  and 
its  prediction,  especially  for  localized 
corrosion.  Extreme  value  statistics  is 
quite  useful  for  assessing  the  maximum 
pit  depth  and/or  the  minimum  time  for 
crack  generation.  The  maximum  pit 
depth  depending  on  the  surface  area 
can  be  evaluated  by  using  the  Gumbel 
distribution  with  the  concept  of  return 
period.  A standardized  procedure  is 
proposed  for  estimating  the  maximum 
pit  depth.  Examples  of  corrosion  failure 


analysis  using  extreme  value  statistics, 
which  were  reported  mainly  in  Japan, 
are  shown.  Accumulated  experiences 
suggest  that  an  appropriate  classifica- 
tion of  data  based  on  the  corrosion 
mechanism  is  required  before  applying 
extreme  value  analysis. 
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1.  Introduction 

Development  of  extreme  value  statistics  and  its 
application  to  various  fields,  including  corrosion, 
has  been  described  by  Gumbel  in  Ref.  [1]  and  [2], 
Evans  is  one  of  the  pioneers  of  modern  corrosion 
science,  and  first  established  the  concept  of  corro- 
sion probability  [3,  4].  Eldredge  [5]  used  extreme 
value  statistics  to  obtain  the  maximum  value  of  pit 
depth  on  an  oil  well  tube  wall  as  a function  of  tube 
surface  area.  Scott  [6]  found  a logarithmic  depen- 
dence of  the  maximum  pit  depth  on  surface  area, 
and  explained  that  dependence  by  referring  to  Trip- 
pet  [7].  Aziz  [8]  and  Eledredge  [5]  discussed  almost 
all  important  points  to  be  considered  for  the  analy- 
sis of  corrosion  pit  data  and  made  use  of  the  return 
period  concept.  This  concept,  originally  introduced 
in  the  fields  of  hydrology  or  meteorology,  is  now 
used  to  obtain  a size  factor  which  makes  it  possible 
to  estimate  the  maximum  pit  depth  in  a large  sur- 
face area  based  on  the  distribution  of  a small  num- 
ber of  pit  depth  data  from  the  small  surface  area 


samples.  In  Japan,  early  review  papers  on  corrosion 
probability  and  extreme  value  statistics  by  Masuko 
[9]  and  Shibata  [10,  11]  contributed  to  the  study  of 
the  extreme  value  statistics  as  applied  to  corrosion 
problems.  Ishikawa  [12,13,14]  and  Imagawa 
[15,16,17]  applied  extreme  value  statistics  to  ana- 
lyze engineering  data.  Kase  [18,  19]  reviewed 
Lieblein’s  paper  [20],  introducing  MVLUE  (mini- 
mum variance  linear  unbiased  estimator)  method 
for  estimating  the  distribution  parameters.  Lieblein 
had  given  the  coefficient  of  MVLUE  up  to  N = 6. 
Recently,  Tsuge  [21]  had  calculated  the  coefficients 
up  to  N = 45,  and  confirmed  that  the  parameters  es- 
timated by  the  MVLUE  method  are  unbiased  and 
efficient  and  are  consistent  with  values  estimated  by 
the  method  of  moments  or  maximum  likelihood 
when  the  sample  size  exceeds  more  than  20.  The 
committee  of  Japan  Society  of  Corrosion  Engineer- 
ing (JSCE)  proposed  a standard  procedure  [22]  to 
estimate  the  maximum  pit  depth  from  the  small 
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sample  size  data  by  using  the  MVLUE  method,  and 
a computer  program  named  EVAN  [23]  was  devel- 
oped. Recently  Laycock  et  al.  [24]  reported  that  a 
generalized  extreme  value  distribution  is  more  con- 
venient for  corrosion  depth  analysis,  because  no 
preliminary  assumption  on  the  type  of  distribution 
is  needed.  An  introductory  book  [25]  by  Kowaka  et 
al.  helped  to  differentiate  extreme  value  analysis 
among  corrosion  workers  in  Japan.  General  back- 
ground on  extreme  value  statistics  is  provided  by 
Ang  and  Tang  [26],  and  Kinnison  [27]. 

2.  Application  of  the  Extreme  Value 
Analysis  to  Corrosion 

In  the  early  1980s,  meetings  and  symposia 
[28,29,30,31]  were  held  in  Japan  for  discussing  the 
basic  principles  of  extreme  value  statistics  as  well  as 
difficulties  and  problems  in  their  application  to  cor- 
rosion. In  Table  1,  several  topics  for  which  the 
Gumbel  distribution  is  applied  are  listed.  Table  2 
includes  cases  analyzed  using  the  Weibull  distribu- 
tion, including  the  exponential  distribution.  Before 
discussing  case  histories,  the  standard  procedure 
[22]  proposed  by  the  committee  is  briefly  explained; 
details  are  available  elsewhere  [22,  32]. 

Table  1.  Examples  of  the  extreme  value  analysis  for  corrosion 
using  the  Gumbel  distribution 


Example  Ref. 


1.  Life  prediction  of  super  heater  tubes  of  the  power  [49] 
plant 

2.  Application  of  the  extreme  value  analysis  to  [37] 

heating  tubes  of  the  boiler 

3.  Estimation  of  the  maximum  amount  of  impurity  [42] 
segregation  in  steel 

4.  Failure  life  estimation  of  SCC  for  Ni  base  alloys  [41] 

5.  Extreme  value  analysis  of  the  corrosion  depth  of  [35] 

the  oil  tank  plate  [36] 

[50] 

6.  Life  prediction  of  heat  exchanger  tubes  [51] 

7.  Eddy  current  examination  system  for  heat  exchanger  [43] 
tubes  with  the  extreme  value  analysis 

8.  Extreme  value  analysis  of  pitting  corrosion  of  heat  [52] 
exchanger  tubes 

9.  Methods  for  the  parameter  estimation  of  the  pit  [16] 
distribution  in  plants 

10.  Ultrasonic  method  for  heat  exchanger  tubes  with  [44] 
the  extreme  value  analysis 

11.  Maintenance  system  for  coated  heat  exchanger  tubes  [45] 

12.  Corrosion  of  steels  in  sea  water  [53] 

13.  Analysis  of  perforation  of  zinc  plating  steels  by  [54] 
extreme  value  statistics 

14.  Fatigue  crack  behavior  of  high  strength  steel  in  [55] 
artificial  sea  water 


Table  2.  Examples  of  the  extreme  value  analysis  for  corrosion 
using  the  Weibull  (exponential)  distribution 


Example  Ref. 


1.  Failure  life  analysis  of  stress  corrosion  [56] 

cracking  of  stainless  steel  heat  exchanger  tubes 

2.  Failure  life  distribution  of  stainless  steels  [57] 

in  high  temperature  and  high  pressure  water 

3.  Effect  of  CaCL  concentration  on  SCC  life  time  [58] 

distributions  of  stainless  steels 

4.  Evaluation  of  SCC  failure  life  of  stainless  steel  [59] 

in  high  temperature  water 


2.1  The  Gumbel  Distribution 

The  procedure  is  proposed  mainly  for  analyzing 
pit  depth  distribution  by  using  the  Gumbel  distribu- 
tion and  the  return  period  in  order  to  estimate  the 
maximum  depth  of  the  larger  surface  area  from 
which  small  area  specimens  are  sampled.  The  Gum- 
bel distribution  is  expressed  as 

F (x)  = exp  (-exp  (-  (x  - \)/a),  (1) 

where  F(x)  is  the  cumulative  probability  of  pit 
depth,  x,  and  A and  a are  the  location  and  scale 
parameters.  The  reduced  variate,  y, 

y=(x-\)/a  (2) 

is  introduced,  and  then 

y = -ln(-ln(E(y)))  (3) 

is  used  for  constructing  the  Gumbel  probability  pa- 
per. Plotting  position  for  the  cumulative  probability 
can  be  calculated  simply  by 

F(y)=  1 — i/(l  +N),  (4) 

where  i is  the  / th  of  the  ordered  value,  x,  in  de- 
scending order  and  N is  the  total  number  of  sample. 
Plotting  y as  a function  of  x yields  a best-fitting 
straight  line;  its  slope  provides  1 /a  and  its  intercept 
(aty  =6)  yields  A.  Instead  of  this  graphical  estima- 
tion of  the  parameters,  more  reliable  estimates  of  a 
and  A can  be  obtained  by  using  the  MVLUE  (min- 
imum variance  unbiased  estimator)  method,  the 
maximum  likelihood  and  the  method  of  moments. 
Among  them  the  MVLUE  method  which  is  dis- 
cussed by  Lieblein  [20]  is  found  to  be  more  efficient 
and  unbiased  for  small  size  samples.  The  MVLUE 
estimator  can  be  calculated  by 
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A =Xflj  (N,  n)xi 

a = Xb,  (TV,  n)x,,  (5) 

where  a,  (N,  n ) and  6,  (N,  n ) are  weights  for  each 
sample  depending  on  the  sample  size,  N,  and  trun- 
cated number,  n,  which  are  tabulated  in  the  table 
given  by  Tsuge  [21,  23]  up  to  TV  = 45.  The  weights, 
A,  B,  C,  of  variance,  V, 

V = a2  (A  ( N,n)y2  + B(N,n)y+C(N,n ))  (6) 

are  also  found  in  the  table  given  by  Tsuge  [21,  23]. 
For  the  pit  depth  distribution,  the  return  period,  T, 
is  defined  as 

T=S/s,  (7) 

where  S is  the  surface  object  (e.g.,  a tank  plate)  to 
be  examined  and  s is  the  area  of  the  small  speci- 
mens which  are  sampled  randomly  from  the  objec- 
tive. The  return  period,  T,  is  in  effect  a size  factor. 
The  mode,  A,  of  the  pit  depth  distribution  for  the 
small  specimen  is  simply  obtained  by  the  MVLUE 
estimators  mentioned  above,  and  the  mode  for  the 
T times  larger  surface,  xmax,  is  given  by 

*max=  A + aln  (T).  (8) 

The  perforation  probability,  P,  of  the  maximum  pit 
through  the  wall  thickness,  d,  is  given  by 

P = 1 — exp  ( — exp  ( — (d  - (A  + aln  (T)/a)).  (9) 

Finally,  the  procedure  [22]  requires  reporting  the 
surface  area  of  the  object,  S,  with  the  small  sample 
area,  s,  providing  the  return  period,  T (-S/s),  and 
the  number  of  samples,  N,  with  data  number,  n , ac- 
tually obtained.  In  addition,  the  original  thickness 
of  the  plate,  d,  and  the  perforation  probability,  P, 
if  needed  are  to  be  stated.  The  above  procedure 
does  not  request  to  check  a goodness  of  fit  of  the 
distribution  obtained  to  the  Gumbel  distribution, 
but  recommends  to  examine  the  fitness  by  the  Kol- 
mogorov-Smirnov  or  chi-square  test  if  needed. 

2.2  The  Weibull  Distribution 

The  third  type  for  the  smallest  value  called  the 
Weibull  distribution 

F (t ) = 1 — exp  ( — ((t  — y)/T|)m)  (10) 

can  be  fitted  to  the  failure  life  distribution  of  stress 
corrosion  cracking  [33,  34]  as  shown  in  Table  2, 


where  y,  77  and  m are  the  location,  scale  and  shape 
parameter,  respectively.  This  third  type  of  asymp- 
totic distribution  for  the  smallest  value  can  be  trans- 
formed to  the  first  type  for  the  largest  value,  that  is, 
Eq.  (1),  by  changing  l-F(t)  to  F(z)  and  by  intro- 
ducing the  following  reduced  variate 

Z = ln(/-y),  z = (X-k)/a.  (11) 

The  same  MVLUE  method  used  for  Eq.  (1)  can  be 
utilized  for  parameter  estimation,  because  the  fol- 
lowing relations  exist  between  the  parameters  of 
both  distributions; 

A = ln(r7 ),  a - Mm  . (12) 

The  above  unified  procedure  for  estimating 
parameters  of  the  Gumbel  and  Weibull  distribution 
was  coded  in  the  computer  program  EVAN  [23]. 


3.  Examples 

Several  examples  are  provided  to  demonstrate 
the  usefulness  of  extreme  value  statistics  for  analyz- 
ing corrosion  problems. 

3.1  Maximum  Pit  Depth  of  Oil  Tank  Plate 

Through  the  1960s  and  1970s  a number  of  oil 
tanks  were  built  in  Japan.  In  the  late  1970s  there 
occurred  frequent  oil  leakages  from  tanks  due  to 
corrosion  failure.  Oil  refinery  or  petrochemical  in- 
dustries were  located  along  the  seacoast  and  oil 
leakage  caused  serious  environmental  damages.  In 
1976,  the  fire  service  law  was  revised  to  enforce  in- 
spection of  the  thickness  of  the  base  and  annular 
plates  of  oil  tanks  every  time  oil  was  evacuated.  On 
these  occasions  extreme  value  analysis  was  applied 
and  found  to  be  a powerful  tool  for  estimating  the 
maximum  pit  depth.  It  is  emphasized  that  data  for 
the  base  plate  and  the  annular  plate  should  be  con- 
sidered separately  because  they  are  characterized 
by  different  corrosion  damage  and  mechanisms. 

The  law  requests  that  plate  thickness  has  to  be 
measured  at  the  corners  of  every  10  cm  square  on 
the  whole  surface  of  the  plate.  This  inspection  pro- 
cedure contributed  greatly  to  reducing  corrosion 
leakage,  but  was  time-consuming  and  costly.  The 
extreme  value  analysis  was  then  studied  intensively 
in  this  field  [35]  [36].  Pit  depth  distribution  sampled 
from  the  whole  base  plate  was  found  to  obey  the 
Poisson  distribution. 
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Araki  et  al.  [36]  found  that  the  largest  value  from 
each  small  square  (5  = 1 m2)  being  randomly  sam- 
pled obeys  the  Gumbel  distribution  as  shown  in 
Fig.  1 [36].  The  slope  and  intersect  of  the  line 
(a  = 0.694  and  A = 1.41)  were  estimated  by  the 
MVLUE  estimators  of  Eq.  (5).  In  this  case,  the  sur- 
face area  of  the  base  plate,  S,  was  1535  m2  and  the 
return  period  or  size  factor  was  calculated  to  be 
T = S/s  — 1535.  The  maximum  depth,  xmax,  was  cal- 
culated by  Eq.  (8): 

xm-M  - 1.41  + 0.694  x In  (1535)  = 6.50  mm 

which  is  shown  also  in  Fig.  1.  These  data  were  ob- 
tained for  a base  plate  which  was  exposed  for  12.6 
years.  Data  for  both  annular  and  base  plate  exposed 
for  7.7  years  were  plotted  in  Fig.  2 [36],  from  which 
the  first  leak  due  to  the  maximum  pit  is  to  be  ex- 
pected after  17.6  years  for  the  base  plate  and  23.5 
years  for  the  annular  plate,  respectively.  The  effect 
of  N,  s,  on  the  estimates  was  examined  and  it  was 
concluded  that  the  MVLUE  method  is  optimal  for 
N < 20  and  the  maximum  likelihood  method  is  reli- 
able for  N > 20. 

3.2  Rupture  of  Heat  Exchanger  Tubes  of  the 

Boiler 

Super  heater  and  economizer  tubes  of  boilers  are 
exposed  to  high  temperature  gases  with  salt  de- 
posits which  cause  severe  corrosion  attack.  Corro- 
sion attack  is  not  uniform,  but  localized  at  several 
sites,  wall  thinning  at  the  localized  site  resulting  in 
burst.  Regular  inspection  is  needed  to  predict  time 
for  replacement  of  the  tube  before  burst.  All  tubes 
have  to  be  examined  for  predicting  exact  time  with 
high  confidence,  but  cost  of  inspection  being  high, 
that  Fukuda  et  al.  [37]  introduced  the  use  of 
extreme  value  analysis  to  supplement  the  inspection 
of  a small  number  of  tubes.  In  Fig.  3,  the  largest  val- 
ues of  wall  thinning  observed  for  14  tubes  are  plot- 
ted on  Gumbel  probability  paper.  The  distribution 
of  wall  thinning  at  every  inspection  time  is  seen  to 
obey  the  Gumbel  distribution  and  the  maximum 
thickness  determined  by  the  return  period  (40 
tubes)  increases  with  operation  time  as  shown  in 
Fig.  4.  A criterion  for  a proper  replacement  time 
has  been  proposed,  which  requires  replacement 
when  the  wall  thickness  reaches  half  of  the  design 
thickness,  tsy.  Risk  of  burst  could  be  avoided  by 
estimating  the  depth  and  noting  the  proposed  crite- 
rion. 


Fig.  1.  The  Gumbel  plot  of  the  maximum  pits  on  the  bottom 
plate  of  the  oil  tank. 


Fig.  2.  Estimated  depth  of  the  maximum  pit  on  the  whole  sur- 
face as  a function  of  operation  years,  and  prediction  of  failure 
life  for  the  bottom  and  annular  plate. 


3.3  The  Pit  Depth  Distribution  of  Steel  Piles  in 
Sea  Water 

Since  the  1970s,  steel  pipes  and  piles  have  been 
used  extensively  in  Japan  for  harbor  construction, 
because  lead  time  for  construction  could  be 
reduced  compared  with  using  concrete.  Recently, 
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Fig.  3.  The  Gumbel  plots  of  the  maximum  thickness  loss  of 
boiler  tubes  used  for  different  operation  times. 


0 1.0  2.0 


Operation  time  (x  104h) 

Fig.  4.  Thickness  loss  of  boiler  tubes  as  a function  of  operation 
time  and  estimation  of  rupture  time. 


corrosion  of  steel  pipes  and  piles  was  found  to  cause 
the  collapse  of  harbor  structures.  Then  corrosion 
damage  of  steel  structures  exposed  in  sea  water  has 
been  inspected  and  analyzed  by  using  extreme  value 
analysis.  Itoh  et  al.  [38]  reported  that  three  differ- 
ent types  of  depth  distribution  were  found  for  steel 
piles  and  plates  depending  on  exposure  time  and 
exposure  location  such  as  water  level  and  deep  sea 
(Fig.  5).  The  type  A distribution  which  exhibits  a 
nearly  straight  line,  was  found  for  uniform  corro- 
sion loss,  its  mean  value  being  below  1.0-1. 2 mm 
thickness.  The  type  C distribution  obeying  the 
Gumbel  distribution  was  observed  for  heavily  local- 
ized specimens.  The  type  B distribution  is  a mixed 
type  of  A and  C distributions.  The  estimated  depth 
using  the  return  period  was  consistent  with  observa- 
tions. 


Thickness  loss  (mm) 

Fig.  5.  Various  types  of  the  distribution  observed  for  steel  piles 
and  pipes  exposed  in  sea  water. 


3.4  Classification  of  Data  Based  on  Corrosion 
Knowledge 

In  any  of  the  cases  mentioned  earlier,  measured 
sets  of  data  is  fitted  by  two  or  three  distributions 
and  must  be  separated  from  each  other  before  the 
analysis  in  order  to  obtain  the  maximum  value. 
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Imagawa  et  al.  [15,16,17,39]  presented  many  cases 
which  require  classification  of  data.  For  example, 
data  for  the  heat  exchanger  tubes  had  to  be  classi- 
fied into  the  inlet  and  outlet  side  samples  because 
corrosion  form  and  its  degree  of  damage  were  dif- 
ferent at  the  two  locations  owing  to  exposure  to  dif- 
ferent temperatures.  For  the  oil  tank,  Imagawa 
observed  that  more  deep  pits  were  formed  on  the 
welding  line  compared  with  other  parts.  He  ob- 
tained the  different  estimated  value  of  the  pit  depth 
for  each  classified  sample.  At  the  present  time,  the 
classification  was  done  on  corrosion  knowledge  and 
experience,  but  it  is  required  to  establish  a proce- 
dure based  on  a common  criterion. 

3.5  Crack  Depth  Distribution  of  Stress  Corrosion 
Cracking 

Stress  corrosion  cracking  is  one  of  the  most  dan- 
gerous corrosion  failure  and  shows  random  occur- 


rence which  is  a very  specific  and  common  feature 
of  materials  fracture.  The  Weibull  distribution  has 
been  known  to  be  quite  useful  to  analyze  the  distri- 
bution of  fracture  strength  of  various  materials  [40] 
and  also  has  been  found  to  be  applicable  for  analyz- 
ing failure  life  distribution  due  to  stress  corrosion 
cracking  [33,  34], 

An  interesting  application  of  the  Gumbel  distri- 
bution for  analyzing  the  crack  depth  distribution 
has  been  reported  by  Tsuge  [41].  The  laboratory  ex- 
periment for  evaluating  the  susceptibility  of  stress 
corrosion  cracking  of  Type  304  stainless  steel  was 
done  by  using  a bent  specimen  of  u-shape.  Bending 
gives  stress  to  the  specimen  and  the  environmental 
condition  of  high  pressure  water  causes  many 
cracks,  which  can  be  revealed  by  sectioning  the 
specimen  after  the  test  as  shown  in  Fig.  6.  Distribu- 
tion of  the  crack  depth  plotted  in  the  Gumbel  prob- 
ability paper  showed  two  lines  with  an  inflection 


Fig.  6.  The  distribution  of  intergranular  corrosion  attacks  and  cracks  observed  for  sensitized 
type  304  stainless  steel  exposed  to  the  BWR  simulated  water  (DO  8 ppm,  250  °C). 
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point  as  can  be  seen  in  Fig.  7.  This  inflection  point 
was  found  to  correspond  just  to  a depth  for  initiat- 
ing the  intergranular  crack.  Thus  the  initiation  of 
the  intergranular  crack  growth  could  be  separated 
from  the  initial  process  of  purely  chemical  inter- 
granular corrosion. 

3.6  Estimation  of  the  Maximum  Segregation  of 
Impurities  in  Steel 

Continuous  casting  of  steel  is  one  of  the  innova- 
tive technologies  achieved  by  the  steel  industry. 
Segregation  and  its  band  which  are  formed  during 
solidification  at  the  center  of  slab  remain  after 
rolling  and  work  as  initiation  sites  for  fracture  phe- 
nomena such  as  lamellar  tear  and  hydrogen  induced 
cracking  (HIC).  The  maximum  amount  of  segrega- 
tion was  found  to  be  related  to  the  above  fracture 


phenomena,  so  that  extreme  value  analysis  was  ap- 
plied for  estimating  the  maximum  amount  of  segre- 
gation in  steel  plate  from  small  area  samples  [42], 
the  concentration  of  impurities  being  measured  by 
using  EPMA.  The  maximum  amount  of  segregation 
thus  determined  can  be  used  to  predict  the  suscep- 
tibility to  lamellar  tear.  It  should  be  emphasized 
that  the  statistical  procedure  for  the  chemical  anal- 
ysis is  mainly  concerned  with  the  mean  and  stan- 
dard deviation  which  assesses  the  reliability  of  the 
measurement,  but  not  with  extreme  values.  In  re- 
cent years,  highly  sensitive  analytical  methods  have 
been  developed,  but  it  is  not  clear  how  to  correlate 
the  data  of  small  area  samples  to  that  of  the  total  or 
bulk  specimen.  The  ratio  of  the  analytical  area  to 
the  bulk  specimen  reaches  almost  to  10_<;,  and  ex- 
treme value  statistics  is  expected  to  be  useful. 


y F 


Fig.  7.  Distributions  changing  with  exposure  time,  the  initial  distribution  corre- 
sponding to  intergranular  corrosion  and  the  second  to  intergranular  cracking. 
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3.7  Non-Destructive  Methods  With  the  Extreme 
Value  Analysis 

Various  types  of  nondestructive  methods  are 
used  for  inspecting  and  examining  corrosion  dam- 
age in  order  to  prevent  failure.  High  sensitivity  and 
resolution  in  time  or  in  space  are  required  for  the 
measurement.  In  addition,  a computer-aided  opera- 
tion becomes  popular,  because  huge  amounts  of 
data  must  be  evaluated.  For  the  heat  exchanger,  a 
thousand  tubes  must  be  checked  and  the  number  of 
data  easily  exceeds  105.  An  eddy  current  sensor  [43] 
and  ultrasonic  sensor  probe  [44]  to  steel  tubes,  and 
an  impedance  sensor  probe  [45]  for  coating  tubes 
have  been  developed  with  the  data  logger  and  the 
extreme  value  analysis  software. 

4.  Discussion 

The  size  effect  on  the  maximum  pit  depth  is 
found  to  be  estimated  with  confidence  by  introduc- 
ing the  concept  of  the  return  period.  Theoretical 
bases  of  the  procedure  have  been  provided  by  ex- 
treme value  theory  [1],  Our  experience  shows  that 
the  pit  depth  distribution  obeys  the  normal  or  expo- 
nential distribution,  which  belong  to  the  exponen- 
tial distribution  family.  Thus  the  maximum  values  of 
pit  depth  extracted  from  the  exponential  family  dis- 
tribution may  reasonably  be  expected  to  obey  the 
Gumbel  distribution.  Thus  the  size  effect  could  be 
rationally  predicted  by  using  the  concept  of  return 
period. 

Evans  [46]  pointed  out,  however,  that  some  cases 
as  observed  by  Wormwell  et  al.  [47]  does  not  obey 
a normal  or  exponential  type  distribution,  but  that 
the  tail  of  the  distribution  is  limited  at  a certain 
depth.  Evans  emphasized  that  such  a limit  is  rea- 
sonable for  the  case  of  anodic  reaction  control  situ- 
ation and  this  limited  depth  gives  a rough  indication 
of  the  greatest  pit  depth  to  be  expected  on  a much 
larger  area.  Evans,  however,  did  not  check  another 
possibility  using  the  type  III  distribution  which  has 
an  upper  limit.  Recently  Laycock  et  al.  [24]  dis- 
cussed that  usefulness  of  the  generalized  extreme 
value  (GEV)  distribution: 

F(x)  = exp  (-(1  -k(x  -u)/a)uk  kx  < =a+uk, 

(13) 

because  the  distribution  subsumes  all  three  types 
with  the  sign  of  a shape  parameter,  k.  When  k is 


zero,  negative  or  positive,  the  distribution  changes 
to  type  I,  type  II,  and  type  III,  respectively.  They 
found  that  the  pit  distribution  on  stainless  steels  in 
acidified  chloride  solution  fits  the  GEV  distribution 
with  k positive,  indicating  that  the  type  III  for  the 
largest  value  could  be  fitted.  The  type  III  distribu- 
tion has  a bound  or  a limit  with  increasing  area,  as 
suggested  by  Evans. 

What  sample  size,  or  what  size  of  specimen  area 
should  be  used  are  questions  from  non-specialists  in 
statistics.  For  this  question,  we  proposed  a proce- 
dure or  criterion  for  choosing  5,  N and  T based  on 
the  variance  given  by  Eq.  (6).  The  surface  area,  S , 
of  the  object  is  given,  and  the  sampling  area,  s,  is 
selected  so  as  to  include  at  least  one  pit.  Then 
T(-S/s)  is  obtained.  Accumulated  data  of  the 
parameters  of  a and  A suggest  [32]  that  the  ratio  of 
a/ A for  localized  corrosion  is  below,  or  not  much 
larger  than,  0.3.  Kinnison  [48]  states  that  the  asymp- 
totic theory  predicts  a constant  ratio  of  0.313  for  all 
extreme  value  distributions.  Then  it  can  be  assumed 
that  the  ratio  a/ A,  is  0.3.  If  we  wish  to  control  vari- 
ance within  (A/3)2,  the  following  relation  can  be  de- 
duced from  Eq.  (6) 

(\/a3)2  =A  (N,n)y2  + B (N,n)y  + C . (14) 

Equation  (14)  can  be  solved  fory  or  T as  a function 
of  N and  al  A,  as  plotted  in  Fig.  8.  When  the  ratio 
of  a/A  can  be  equated  to  0.3  as  discussed  before,  a 
suitable  number  of  samples  can  be  found  for  a given 
return  period,  T.  From  this  figure,  the  required  size 
of  samples  is  A = 30  for  T = 1200,  or  N - 20  for 
T = 274  and  so  on.  This  figure  is  approximately  the 
same  as  what  was  observed  empirically. 

5.  Conclusions 

Extreme  value  statistics  has  been  found  to  be  a 
powerful  tool  for  estimating  the  maximum  value  of 
localized  corrosion  depending  on  the  surface  area. 
Accumulation  of  data  and  experience,  however,  re- 
veals that  statistics  is  less  important  than  corrosion 
experience  and  knowledge  for  obtaining  a reason- 
able estimation;  measured  data  must  be  classified 
based  on  the  form  of  corrosion  damage  and  its  de- 
gree before  the  analysis.  Properly  classified  data  is 
found  to  provide  a very  reasonable  value.  Nonde- 
structive methods  for  measuring  wall  thickness  with 
various  types  of  sensors,  combined  with  extreme 
value  analysis,  have  been  developed  in  recent  years. 
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Fig.  8.  Calculated  curves  of  the  optimum  condition  to  choose  T 
and  N at  various  value  of  a! A. 
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